Quantum filtering one bit at a time 
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In this paper we consider the purification of a quantum state using the information obtained from 
a continuous measurement record, where the classical measurement record is digitized to a single 
bit per measurement after the measurements have been made. Analysis indicates that efficient and 
reliable state purification is achievable for one- and two-qubit systems. We also consider quantum 
feedback control based on the discrete one-bit measurement sequences. 
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In classical systems, filtering a sequence of measure- 
ments in order to reconstruct the behavior of a classical 
object is well established The measurements repre- 
sent information that has been obtained about the cur- 
rent state and dynamical evolution of the system. The 
state constructed using this data represents the 'best 
guess' of the state of the system given the measurement 
record. This state-estimation process is often referred to 
as tracking In quantum systems, an equivalent ap- 
proach has been developed using continuous weak mea- 
surements 0, This approach has shown a range of 
interesting results using the reconstructed (or partially 
reconstructed) quantum state to apply controls to mod- 
ify the subsequent quantum evolution Q. Of particular 
interest here is rapid state purification introduced by Ja- 
cobs [|| and later studied by other authors [fHH . In rapid 
state purification, the rate at which a state can be puri- 
fied can be modified by applying simple controls. 

A major difference between classical state-estimation 
filters and quantum filtering is that, in quantum systems, 
the act of measurement has an irreducible effect on the 
system being measured. A projective (von Neumann- 
type) measurement is an extreme example, leaving the 
quantum system in an eigenstate of the measured opera- 
tor. In continuous weak measurement models, the quan- 
tum system of interest is considered to be coupled to an 
environment. Projective measurements are taken on the 
environment, which contains information on the system 
of interest through the coupling. By averaging over the 
environmental degrees of freedom, it is possible to form 
an equation for the evolution of the quantum system of 
interest, which contains Hamiltonian evolution that has 
been perturbed by the effect of the environmental mea- 
surements. The average evolution of the quantum system 
in the presence of the environment is given by the master 
equation, and the particular evolution for a given set of 
measurements (a realization) is represented by a stochas- 
tic master equation (SME) Q . The conditioned evolution 
for a given measurement record is often referred to as an 
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'unraveling' of the master equation Q. Such unravelings 
may be continuous or discontinuous ('quantum jumps'), 
depending on the nature of the interaction and the mea- 
surement process. They represent quantum filtering pro- 
cesses similar to a classical tracking filter. The measure- 
ment record generated by a continuous weak measure- 
ment is a classical analog signal. The signal is typically 
very noisy (quantum noise generated by the back action 
from the measurement interaction) and integration over 
discrete time steps allows the SME to be calculated nu- 
merically to reconstruct the quantum state. 

This paper is motivated by two factors. Firstly, it 
has been shown recently that for a discontinuous/jump 
unravelling of the master equation, there is a minimum 
number of classical states that are required to track the 
state of a quantum system with D-energy levels In 
the case of a qubit, it was demonstrated that a single 
classical bit (i.e. two classical states) was sufficient to 
track the qubit's state - that is, if one knew the classical 
state at any point in time, one would be able to infer 
the state of the qubit as it jumps between states. Sec- 
ondly, it is known in signal processing that if a classical 
signal is sent through a number of parallel (thresholded) 
channels, each of which is corrupted by Gaussian noise, 
then the optimum number of thresholds that can be used 
to discretize this signal reduces to one as the noise level 
increases (i.e. one bit per channel) [l(|. This is related 
to a classical nonlinear effect: supra-threshold stochastic 
resonance (llj . 

In the cases considered here, parallel channels (i.e. 
multiple measurement interactions) would introduce ad- 
ditional noise through the measurement back action. In- 
stead, we use the fact that the underlying signal is slowly 
varying with respect to the measurement/purification 
rate to perform an (implicit) temporal average over the 
noise. The time required to complete the purification 
is dependent on the coupling to the environment (the 
measurement strength), which is far longer than the in- 
cremental measurements (or samples) used in the SME. 
We show that a single threshold (a one-bit record, OBR) 
in the SME is sufficient to purify the state of one or two 
qubits, with or without feedback control based on local 
unitary operations. Crucially, the thresholding of the sig- 
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nal is performed after the classical measurement record 
is generated so it does not add to the back action on the 
quantum system. We demonstrate that the fidelity of 
one-bit filtering can be very good, reproducing classical 
and quantum correlations between qubits (as measured 
by the quantum discord [Hj]). The importance of these 
results rests in the potential for experimental realization 
of quantum filtering - reducing a continuous time analog 
signal to a discrete-time single bit measurement sequence 
provides a significant reduction in the data required to 
reconstruct the quantum state of the system of interest. 

In Jacobs' rapid state purification |5(, a qubit starts 
in a completely mixed state (i.e. an unknown state) and 
undergoes controlled rotations to rotate its Bloch vec- 
tor onto the plane orthogonal to the measurement axis. 
These controlled rotations alter the rate of increase of 
the average purity by a factor that approaches two as 
the purity of the state P — > 1 dctcrministically (where 
P{p) = Tr(p 2 )). Alternatively, if feedback is used in the 
opposite manner, to rotate the Bloch vector towards the 
measurement axis, purification is stochastic and the time 
required to reach agiven purity can be reduced relative to 
Jacobs' protocol jg]. (Jacobs considered discrete 'bang- 
bang' controls @ , but this is distinct from the discretized 
measurement and filtering process discussed here). 

In this paper, we use two stochastic master equations: 
one SME to generate a measurement record and another 
to estimate the state of the system conditioned on the 
measurement record. The first SME (SMEl) is initialized 
in a pure state and remains in a pure state (po). If an 
experiment were available, SMEl would be redundant, 
its main role is to generate the measurement record. The 
secondary role of SMEl is to provide a reference against 
which we can compare our estimate of the state. SME2 
generates the estimated state. SME2 is initialized in a 
completely mixed state, which will purify gradually as 
information is extracted. All feedback controls are based 
on SME2 and the estimated state only. 

Continuous weak measurements are not restricted to 
Hcrmitian operators but - for simplicity - we will assume 
that the measurement operator is Hermitian. For Her- 
mitian weak measurements dj = y {2kjK)y,j 



c-, where 



kj is the strength of the jth measurement interaction, 
the stochastic master equation (SMEl) for the density 
matrix, po, is given by |3j, 

dp a = - l - po k i \Vj » \Vj > Po\] dt 

j 

+ Y (yj p ° + p°yj ~ 2 (&)o po) dw 3 w 



where H is the Hamiltonian of the system, dt is an in- 
finitesimal time increment and dWj is a real Wiener in- 
crement (such that dWj = and dWjdWy = Sjjidt). 
The corresponding continuous measurement records arc 



To solve this problem numerically, finite time steps At 
are integrated (Euler method) to find a discrete mea- 
surement Ayj(t). This is a classical analog signal, which 
represents the sampled signal that would be derived from 
an experiment. For the one-bit measurement record, the 
discrete samples are replaced with, 



AYj(t) = y/(hAt)sgn{A yj ) 

where sgn(x) = ±1 depending on the sign of x. The 
estimated state, p c , is constructed using SME2, 

Ap c = -Mf,/9 c At-^kj [y^ly^pcWAt 
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+ Y, V^j (VjPc + Pcilj - 2 (y j ) c Pc) AW,-(3) 



where AWj = AYj(t)/Vfr- y/8k~ {yj) c At. This is anal- 
ogous to the classical innovation, the difference between 
the actual measurement and that predicted using the es- 
timated state [l|. 

Here we will consider two specific cases: one qubit with 
a single measurement along the Z-axis, and two coupled 
qubits with a measurement of ZZ = o~ z g) a z (where o~ r 
is a Pauli matrix with r = x,y,z). In each case, we 
will assume that the qubits have a background Hamilto- 
nian, which corresponds to a rotation about the X-axis 
with an oscillation frequency ojq. Such Hamiltonians are 
common in solid state qubits, where they correspond to 
tunneling interactions. Local controls are used because 
they are seen as a tractable way to preserve entangle- 
ment for quantum processing applications [l3[. We will 
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use the fidelity F = F(p ,p c ) = |Tr [^J ^/p~cPo^fp~c\ \" [li 
to characterize the errors in our estimate p c of the true 
state po- For the two qubit case, we calculate classical 
correlations C(p) [l2j], 



C(p) = S(p B ) 



(2) 



and the quantum discord T>(p) = S(p A ) + S(p B ) — 
S(p)—C(p) [l2|, where S(p) is the von Neumann entropy 
S(p) = —Tr(p log 2 p) and p A and p B are the reduced 
density matrices for the individual qubits, and the mini- 
mum is taken over the set of projectors on qubit A, {n^}. 
S , (/o B |^4jn A }) is the entropy of qubit B under the action 
of the projectors on qubit A, and p a is the probability 
associated with the result obtained for that projector on 
qubit A. Orthogonal projectors are used here because, 
although the minimum value that they give is not always 
the true minimum, the value of the quantum discord is 
very close to the true value [TB ]. 

For the single qubit case we have a Hamiltonian H = 
fiwoc x /2, an initial state p c (0) = 1/2 and a single mea- 
surement operator t\ = \J (2k 1 h)o- z . In Figure 1, we 
show the average fidelity as a function of time for the 
case where no controls are applied and the single-qubit 
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FIG. 1: (Color online) Average fidelity of one qubit as a func- 
tion of time with no controls: full record (red-dash) and OBR 
(blue-solid). Inset: an example run, showing fidelity for one 
measurement record. 
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FIG. 2: (Color online) Average fidelity of one qubit as a 
function of time with controls applied to rotate Bloch vec- 
tor onto measurement axis: full record (red-dash) and OBR 
(blue-solid). Inset: Average fidelity after 50 cycles as a func- 
tion of the angle to the measurement axis (in X-Z plane) for 
OBR. 



Bloch vector rotates freely about the X-axis. The (an- 
gular) frequency luq is fixed and other parameters are 
scaled relative to luq: k\ = 0.005o;o. SME2 is integrated 
numerically using the Euler method with 10000 steps per 
cycle and averaged over 500 measurement records. In the 
limit At — > dt, the first two moments of the OBR and 
the full record are identical (mean zero, variance hdt) 
and it is possible to show that the two records will pu- 
rify at the same rate, on average. The analytical results 
in Q also apply for the OBR: the average time to ap- 
proach unit purity is one-quarter of the time to reach an 



average purity of unity. For finite time steps the OBR 
purifies slightly slower than the full record (approx. 20- 
30% longer). However, it is clear that, although the SME 
using the OBR does estimate the underlying state of the 
system reasonably well, the average fidelity for the state 
generated by the OBR saturates at a fixed level < 1. 
In each run of the SME, the fidelity approaches one but 
suffers from occasional drops in fidelity as the errors in 
the state estimated by the OBR accumulate, see Fig- 
ure I (inset). Approximately five to ten thousand OBR 
sample points are required per qubit cycle, longer time 
steps can make the filter unstable, with purities growing 
significantly above one. The number of sample points re- 
quired for a stable filter is related to the qubit frequency 
and the measurement strength. For larger k values, more 
frequent sample points are needed but purification takes 
less time in total. 

To improve the fidelity of the final state estimate it is 
necessary to apply feedback controls to rotate the Bloch 
vector towards the measurement axis at each time step. 
This is the feedback protocol studied in Q. Figure 2 
shows the effect of adding the feedback control. The fi- 
delity of the final state for the OBR approaches one as 
the system purifies, albeit slightly slower than for the full 
record. This is true whether there is Hamiltonian evolu- 
tion or not. In the case where there is no Hamiltonian 
evolution, the protocol in @ reduces to a classical (no 
feedback) measurement, but in either case, the OBR al- 
lows the state to be purified and the fidelity of the final 
state approaches one in the long time limit. The inset 
in Figure 2 shows the dependence of the final (average) 
fidelity on the angle between the measurement axis and 
the vector towards which the Bloch vector is rotated. 
This shows that the fidelity of the final purified state is 
relatively insensitive to small errors in the feedback con- 
trols - the final fidelity is very close to one for all angles 
< 20°. 

For two qubits, the situation is more complex with 
the need to estimate the combined state and (ideally) 
to retain the ability to generate entangled states. Start- 
ing with a two qubit Hamiltonian, Ha+b = ttidoi&x ® 
I) /2 + Hujo(I <g) cr x )/2, there are two main ways to gener- 
ate an entangled state from a separable initial state: an 
entangling Hamiltonian or an entangling measurement. 
An interaction of the form H int = n(a z <£> a z ) can gen- 
erate entanglement. When combined with two local Z 
measurements and no controls, the full record purifies 
the state effectively but the OBR generates an unsta- 
ble state estimate which occasionally drifts away from 
the true state, as with the example shown in Figure 1. 
Using local controls for the individual (reduced) density 
operators, rotating onto the respective iv-axes using ei- 
ther record will purify the two-qubit state with a fidelity 
that tends to one in the long time limit. However, this 
has the disadvantage of removing entanglement between 
the two qubits, because the resultant state is separable 
with each qubit in an eigenstatc of Z . The alternative 
approach is considered here, where entanglement is gen- 
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FIG. 3: (Color online) (a) Average fidelity of the two qubit 
state as a function of time with single ZZ measurement record 
and no controls: full record (red-dash) and OBR (blue-solid), 
and (b) classical correlations (red-cross) and quantum discord 
(blue-circle) of the two qubit state as a function of time with 
single ZZ measurement record and no controls: full record 
(dot) and OBR (solid). 

erated by the measurement. This is similar to the case in 
Q , where a ZZ measurement was used to purify and en- 
tangle two qubits. In Q, the two qubits were separated 
into two encoded qubits and an analogue of Jacobs' pro- 
tocol based on local operations was used to purify and to 
protect entanglement by manipulating the system into a 
decoherence free subspacc. Here, Jacobs' protocol would 
purify the combined state at the expense of the fidelity 



of the final state - as seen in the case of the single qubit 
above - and the alternative approach based on [6| is not 
available with purely local operations. Instead, we opt 
to use no feedback in this case, relying on the inherent 
AT-rotations and a single (entangling) ZZ-measurement. 
Taking an initial state, p c (0) = (I <8> -0/4, k\ ~ 0.005o;o, 
and using 20000 steps per cycle (the same data rate per 
qubit as in the previous example - and subject to the 
same requirement on data rate for a stable filter), the 
purity of the estimated state approaches one for both 
the full record and the OBR. The OBR purifies the state 
slower than the full record, but the fidelity does tend to 
one in the long time limit. The OBR also captures the 
classical and quantum correlations very well (see Figure 
3(b)). Both the full record and the OBR reproduce the 
classical correlations very rapidly, but the quantum cor- 
relations indicated by the quantum discord, D(p), are 
more difficult to estimate using the OBR. 

In this paper, we have demonstrated that quantum 
filtering can be performed using a one bit measurement 
record. Single qubit states and two qubit entangled states 
undergoing continuous weak measurement can be esti- 
mated accurately from a single bit per time step. This 
represents a significant reduction in the data required to 
use quantum filtering methods with only minor reduc- 
tions in filter performance. The one-bit measurement 
record can provide an accurate state estimate for one 
qubit when local feedback is applied, and it has been 
shown that the estimates are insensitive to small feed- 
back errors. For two qubits, it is possible to purify the 
combined state with local measurements and local feed- 
back, but this is at the expense of entangled states. If 
an entangling measurement is used, it is possible to use 
the one-bit record without feedback to generate entan- 
gled states with fidelities approaching one. The one-bit 
record also provides accurate estimates of the classical 
and quantum correlations. 
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